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Abstract 

The reconstruction efficiency of photons and neutral pions is measured using the 
relative yields of reconstructed B+ —)> J/rj; K*+ (—)> and B+ —)> J/ij; K+ decays. 

The efficiency is studied using the data set, corresponding to an integrated lumi¬ 
nosity of 3 fb~^, collected by the LHCb experiment in proton-proton collisions at 
the centre-of-mass energies of 7 and 8 TeV. 


1 Introduction 

The LHCb detector is a single-arm forward spectrometer at LHC covering the 
pseudorapidity range 2 < rj < 5, designed for the study of particles containing b or c 
quarks. Many of the decays under study have photons and neutral pions in the hnal 
state, for example [2]-^. The knowledge of photon reconstruction efficiency in a wide 
range of photon kinematic is a crucial ingredient for such kind of measurement [^, [^. 
This manuscript describes the study of photon reconstruction efficiency using large sam¬ 
ple of B+ —)■ J/ij) K*"*" decays and known ratio of branching fractions of B+ —)■ J/ij) K*"*" and 
B"*" —)■ J/ipK’*' decays [^. A correction factor for the reconstruction efficiency gained from 
simulation, is obtained through the detailed comparison between simulation and data. 

The photon reconstruction efficiency is determined in data using comparison between 
B”^ -A J/i1jK*+ (—;■ K’''7t°) and B"^ -A J/ipK’'' signal yields. The correction factor to 7t° re¬ 
construction efficiency is obtained as the following: 

oo„_ 'fBLj/tK+ 

I." X B (B+ ^ J/4) K-+ K+tiO)) • ^ ' 
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where N is the event yield, is the total efficiency, determined using simulation, and 
the ratio of branching fractions is known to be 


BiB+ 3/^\)K+) 

B (B+ ^ J/\1;K*+ K+7t0)) 



X (1.39 ±0.09) 



Since these decays have similar topology, the large cancellation of various systematic 
uncertanties, after imposing similar selection and trigger criteria is expected. 


2 LHCb detector and data sets 


The LHCb detector includes a high-precision tracking system consisting of a silicon-strip 
vertex detector surrounding the pp interaction region [^, a large-area silicon-strip detector 
located upstream of a dipole magnet with a bending power of about 4 Tm, and three sta¬ 
tions of silicon-strip detectors and straw drift tubes placed downstream of the magnet. 
The tracking system provides a measurement of momentum, p, with a relative uncer¬ 
tainty that varies from 0.5% at low momentum to 1.0% at 200GeV/c. Different types 
of charged hadrons are distinguished using information from two ring-imaging Cherenkov 
detectors 10 . Photon, electron and hadron candidates are identihed by a calorimeter sys¬ 


tem consisting of a scintillating-pad detector (SPD) which is designed to identify charged 
particles, followed by a wall of lead and the preshower detectors (PS), an electromag¬ 
netic calorimeter (ECAL) and a hadronic calorimeter (HCAL). Muons are identihed by a 


system composed of alternating layers of iron and multiwire proportional chambers 11 


The calorimeter system performs several tasks, providing the hrst level trigger with 
high transverse momentum photon, electron and hadron candidates, measuring their ener¬ 
gies and positions and performing the separation between photons, electrons and hadrons. 
All calorimeters use the same energy detection principle: scintillation light is transmitted 
to a photomultiplier tube (PMT) by wavelength-shifting (WLS) hbres. The single hbres 
for the SPD/PS cells are read out using multianode photomultiplier tubes (MAPMT), 
while the hbre bunches in the ECAL and HCAL modules require individual phototubes. 
The SPD/PS, ECAL and HCAL adopt different cell dimensions in the inner and outer 
parts of the calorimeter since the hit density varies by two orders of magnitude over the 
calorimeter surface. A segmentation into three different sections has been chosen for the 
ECAL and projectively for the SPD/PS. Civen the dimensions of the hadronic showers, 
the HCAL is segmented into two zones with larger cell sizes. 

The electromagnetic calorimeter employs ’’shashlik” technology of alternating scintil¬ 
lating tiles and lead plates. The energy resolution of the ECAL modules was determined 
at the test beam. The parametrisation cte/E = a/\/E©b©c/E (E in CeV) is used, where 
a, b and c stand for the stochastic, constant and noise terms respectively. Depending on 
the type of module and test beam conditions the stochastic and constant terms were mea¬ 
sured to be 8.5% < a < 9.5% and b ~ 0.8%. The noise term corresponds to c = 0.003 
CeV. 

Energy deposits in ECAL cells are clusterised applying a 3 x 3 cell pattern around local 
maxima of energy deposition. Consequently the centers of the reconstructed clusters are 
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always separated at least by one cell. If one cell is shared between several reconstrncted 
clnsters, the energy of the cell is redistribnted between the clusters proportionally to the 
total cluster energy. The process is iterative and it converges rapidly due to relatively small 
ratio of the Moliere radius (3.5 cm) to the cell size. After the redistribution of energy 
of shared cells, the energy-weighted cluster moments up to the order 2 are evaluated 
to provide the (hypotheses-independent) cluster parameters, namely, the total energy, 
the energy barycenter position and the two-dimensional energy spread matrix. Neutral 
clusters are identihed as those clusters that do not match to charged tracks. The photon 
energy, Ec, is evaluated from the total cluster energy, £ecal, and the reconstructed energy 
deposit in the PS, £ps, as follows: Ec = a x £ecal + /^ x £ps. The parameters a and 
(3 account for energy leakage. The value of parameter a depends on the position of the 
cluster on the calorimeter surface while f3 is estimated afterwards to give full account of 


energy samples 12,13 


The calibration of the electromagnetic calorimeter is performed through several steps. 
The energy resolution of each ECAL module has been determined at a test beam 14 to 
be in agreement with the designed value of ^ 


= ^© 1 % 


15 


The photoelectron 


multipliers gains are determined with the help of an LED system installed in the calorime¬ 
ter. The hnal step of the calibration procedure is performed with several methods using 
the data coming from the experimental setup. First the energy flow method is applied, 
which is based on the idea that transverse energy flow measured by the ECAL is a lo¬ 
cally smooth function 16 . Further the calibration with neutral pions is performed. This 
method consists in determining the neutral pion peak position in 7t° —>■ yy decays ^18 . 

The performance of high energy photon reconstruction is illustrated by the recon¬ 
structed B° —)■ K*°y and B° —)■ (jjy mass distribution 19 
for this radiative decay, 92 MeV/c^, 


The mass resolution obtained 


is dominated by the ECAL energy resolution. Neu¬ 
tral pions with low transverse energy are mostly reconstructed as a resolved pair of well 
separated photons. A mass resolution of 8 MeV/c^ is obtained for such neutral pions. 

In 2011(2012) LHCb was taking data in proton-proton collisions at a centre-of-mass 
energy of a/s = 7 TeV (8 TeV). The collected data samples correspond to an integrated 
luminosity of 1 fb~^ in 2011 and 2 fb~^ in 2012. The average number of visible pp 
interactions per bunch crossing, p, is 1.5 for both periods, which is four times larger 
than the designed value. The instantaneous luminosity, £, was increased during 2011 up 
to 4 X 10“^^ cm“^c“^ and during 2012 was kept constant at this value which is two times 
larger than the designed value. The excess of /i and C in comparison with the design 


values of 2 X 10“'^^ cm“^c“^ leads to higher event multiplicities, which also implies larger 
number of tracks and higher number of clusters in the electromagnetic calorimeter and, as 
consequence, higher number of overlapping clusters. Due to the difference in conditions of 
the data taking, the correction factors are calculated separately for 2011 and 2012 periods. 
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3 Event selection 


The decays B+ —)■ J/t|) K*+ followed by K*+ —)■ K+7t° and B'*' —)■ J/ij; K+ are reconstructed 
using the J/t|) —>■ decay mode. The signal candidates are selected with the help 

of cuts applied to various kinematic parameters, particle identihcation and decay recon¬ 
struction quality variables. Most selection criteria are common for the decay channels, 
except those which are related to the photons selection. 

The track quality of reconstructed charged particles is ensured by requiring the per 
degree of freedom, y;^j,/ndf, provided by the track £t to be less than three. Well-identihed 
muons are selected by requiring that the difference of logarithms of the muon and hadron 
hypothesis likelihood, provided by the particle identihcation detectors 20 , is 

larger than zero. 

Pairs of oppositely-charged tracks identihed as muons, each having a trans¬ 
verse momentum greater than 550 MeV/c are combined to form J/r|; —)■ 
candidates. The invariant mass of dimuon combination is required to be be¬ 
tween 3050MeV/c^ and 3150MeV/c^ for J/ij) candidates. Only events explicitly triggered 
by the selected J/ip candidates at all trigger levels have been considered in the analysis. 

Photons are reconstructed using the electromagnetic calorimeter 12,13 and identihed 


using a likelihood-based estimator, constructed from variables that rely on calorimeter and 
tracking information. Candidate photon clusters must not be matched to the trajectory 
of a track extrapolated from the tracking system to the cluster position in the calorime¬ 
ter. Further photon quality rehnement is done using information from the PS and SPD 
detectors. The photon transverse momentum is required to be greater than 250 MeV/c. 

For the photon identihcation a variable is used, which combines information from the 
calorimeter and tracking systems in order to reduce the background from hadrons, elec¬ 
trons and merged 7r° (those which daughter photons create largely overlapping clusters in 
the ECAL). The value of this variable, CL, is required to be more than 2%. The invariant 
mass of the photon pair is required to be within 30 MeV/c^ around the nominal mass. 

Charged kaons are selected by requiring the transverse momentum to be in excess 
of 450 MeV/c. Loose kaon identihcation criterion is applied, the corresponding neu¬ 
tral network variable, Pk, is required to be greater than 0.1. The invariant mass of 
the K+7r° combinations is required to be in range ± 75MeV/c^ around the nominal 
K*+ mass [^. 


The B+ candidates are formed from pairs. A reht 


21 


of the B'*' candidate 

is performed taking into account primary vertex pointing constraint and J/ij) and 71° mass 
constraints. The reduced for this procedure, xlt/irdf, is required to be less than hve. 
The transverse momentum of the B+ candidate is required to be larger than 3 GeV/c. To 
reject possible contributions to B+ —)■ J/rpK*"'' decay from B"*" J/tJjK’'' decays with two 
additional random soft photons, the invariant mass of the J/ij) K"*" combination is required 
to be outside a ±25 MeV/c^ mass window around the known B"*" mass. Finally the decay 
time, cr, of B+ candidates is required to be in excess of 200 pm. 
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4 Signal yields 


The invariant mass distributions for the simulated B+ —>■ J/rp K"*" and B’*' —>■ J/r|; K*"*" 
events, which have passed the selection, are shown in Fig. [1} The mass distributions 
for the B+ —>■ J/t|) K’*' and B'*’ —)■ J/i|) K*"*" candidates selected from the data are shown 
in Fig. The diphoton invariant mass distribution in data obtained using the sPlot 
technique is shown in Fig. A clear peak from the t\P —)■ 77 decays is seen in 
the distribution. Its position is found to be 134.9 ± O.lMeV/c^ and the resolution is 


8.3±0.1MeV/c2. 

The signal yields are determined using an extended unbinned maximum likelihood 
£t. The signal component for both decays is modelled by a double-sided Crystal Ball 
function 23 . The tail parameters of the signal function are hxed using simulated events, 
whereas the mean and resolution of the B+ candidate are allowed to vary in the £t. The 
combinatorial background is modelled by an exponential function. For both hnal states 
the htted position of the B"*" peak is consistent with the known B+ mass and the mass 
resolution is in agreement with the values observed in simulated samples. 


5 Efficiencies 


The efficiencies and their ratios are estimated using simulation. For the simulation, pp 
collisions are generated using Pythia 6 and Pythia 8 24 with a specihc LHCb conhg- 
uration 25 . Decays of hadronic particles are described by EvtGen 26 in which hnal 


state radiation is generated using Photos 27 . The interaction of the generated parti¬ 


cles with the detector and its response are implemented using the Geant4 toolkit 
as described in Ref. 


29 . The digitized output is passed through a full simulation of both 


the hardware and software trigger and then reconstructed and selected in the same way 
as the data. 

Ratios of efficiencies between B+ —>■ J/t|)K*+ (—;■ K+7t°) and B+ —>■ J/r|jK''' decay modes 
are calculated separately to obtain the overall efficiency ratio. The following efficiency 
ratio is obtained for 2011 data 


^total 

^B+^J/4)K+ 

ptotal 


13.32 ±0.18, 


and for 2012 data 


^total 

^B+^J/4)K+ 

ptotal 

^B+^J/i[jK*+ 


13.57 ± 0 . 12 , 


where the uncertainty is due to the limited size of simulated samples. 
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Figure 1: Simulated mass distributions for the selected B+ —J/4)K+ (left) and 
B+ —)■ J/i[)K*+ (right) candidates (a,b) for 2011 and (c,d) 2012 data. Red curves rep¬ 
resent the total £t function. Dashed purple and dashed blue curves represent signal and 
background respectively. 


6 Systematic uncertainties 

Since the decay products in each channel have similar kinematic properties, most uncer¬ 
tainties cancel in the ratio, in particular, those related to the muon and J/ij) reconstruction 
and identihcation. The remaining systematic uncertainties are discussed below. 

Systematic uncertainties related to the £t model are estimated using alternative models 
for the description of the invariant mass distributions. The following functions are used 
as alternative: 

• a sum of a Student-t distribution with all parameters hxed from the simulation for 
the signal and exponent for background; 
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• a sum of modified Novosibirsk function with all parameters fixed from the 
simulation for the signal and exponent for background; 

• a sum of Crystall Ball function for the signal and second order polynomial for 
background. 

For each £t model, the ratio of event yields is calculated. The maximum observed devia¬ 
tion from the baseline model in the ratio of yields is taken as systematic uncertainty. The 
obtained uncertainty is 0.2 %. 

The uncertainty due to the limited number of simulated events vary in a range from 
0.8% to 4.5%, depending on the requirements on the photon transverse energy and neutral 




mj/^i,K+ [GeV/c^] 


mj/pK*+ [GeV/c^] 



Figure 2: Invariant mass distributions for the selected B'*' —J/4)K+ (left) and 
B+ —)■ J/i|)K*+ (right) candidates (a,b) for 2011 and (c,d) 2012 data. Red curves rep¬ 
resent the total £t function. Dashed purple and dashed blue curves represent signal and 
background respectively. 


7 



















Figure 3: Invariant mass distributions for the selected 7 r° —77 candidates. 

pion transverse momentum. The corresponding systematic uncertainty is assigned in each 
bin of the photon transverse momentum individually. 

The trigger efficiencies for events with J/t|) —)■ p.+ p~ produced in beauty hadron decays 
is studied and a systematic uncertainty 1.1 % is assigned based on the comparison of the 
ratio of trigger efficiencies for high-yield samples of B"*” —)■ J/i|)K+ and —)■ r[)(2S)K+ de¬ 
cays on data and simulation [^. 

The acceptance efficiency is calculated separately for different magnet polarities and 
for different data taking periods. The obtained individual acceptance efficiencies for the 
B"*" —)■ J/ij) (K*"*" —)■ K’'‘ 7 h’) and B"*" —)■ J/il;K+ decays and their ratios are represented in 
Table The average of four ratios for 2011 (1.1040 ± 0.0015) is chosen as the resulting 
ratio of acceptance efficiencies (uncertainty is statistical only). The maximum deviation 
from the average is taken as systematic uncertainty and equals 1.4 %. The average of four 
ratios for 2012 (1.0998 ± 0.0015) is chosen as the resulting ratio of acceptance efficiencies 
(uncertainty is statistical only). The maximum deviation from the average is taken as 
systematic uncertainty and equals 0.2 %. The observed difference in the efficiency ratios 
is conservatively taken as an estimate of the systematic uncertainty and is 1.4 % for 2011 
data and 0.2 % for 2012 data. 

The agreement between data and simulation has also been cross-checked by varying 
the selection criteria. No unexpectedly large deviation is found and no contribution is 
taken as systematics. 

The contribution to B+ —)■ signal yield from B'*' —J/ipTf'' decay mode, when 

pion is misidentihed as kaon, is found to be negligible (lower than 0 . 001 %). And the 
contribution to B+ —)■ J/4)K*+(—)■ K+ 7 r°) signal yield from B+ —)■ J/t|)P+(—)■ 71+71°) decay 
mode, when pion is misidentihed as kaon, is also hnd to be negligible (around 0 . 02 %). 

The summary of systematic uncertainties is presented in Table The total systematic 
uncertainty is 1 . 8 % for 2011 and 1 . 1 % for 2012 . 
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Table 1: Acceptance efficiencies with their ratios obtained from simnlation. Uncertainties 
on the efficiencies and their ratios are statistical only and reflect the simnlation statistics. 


Period 

Pythia 

Magnet 

B+ ^ J/r|;K+ 

B+ ^ 7/^\>K*+ 

gen&acc / gen&acc 

%+-?• J/^l) K+ / ^B+-5- J/il) K*+ 

2011 

6 

Up 

0.1550 ±0.0003 

0.1396 ±0.0003 

1.111 ±0.003 



Down 

0.1547 ±0.0003 

0.1401 ±0.0003 

1.104 ±0.003 


8 

Up 

0.1636 ±0.0003 

0.1491 ±0.0003 

1.097 ±0.003 



Down 

0.1640 ±0.0003 

0.1485 ± 0.0003 

1.104 ±0.003 

2012 

6 

Up 

0.1578 ±0.0003 

0.1432 ± 0.0002 

1.102 ±0.003 



Down 

0.1573 ±0.0003 

0.1431 ±0.0002 

1.099 ±0.003 


8 

Up 

0.1665 ±0.0005 

0.1518 ±0.0002 

1.097 ±0.003 



Down 

0.1675 ±0.0004 

0.1522 ±0.0003 

1.101 ±0.003 


Table 2: Relative systematic uncertainties for the correction factor of the reconstruction effi¬ 
ciency. The total uncertainty is the quadratic sum of the individual components. 


Source 


Fit model 
Trigger 
Acceptance 
Total 


Uncertainty [%] 

2011 2012 

Vs = 7 TeV Vs = 8 TeV 
0.2 
1.1 

1.4 0.2 

1.8 1.1 


7 Correction factor 

The correction factor is calculated according to Eq. For 2011 data its value is calculated 
to be 

riX = (103.2 ± 2.6 (stat) ± 2.3 (syst) ± %, 

and for 2012 data is calculated to be 

riX = (105.9 ± 1.8 (stat) ± 1.6 (syst) ± %, 

where the third uncertainty is related to an uncertainty of the known branching fractions 
of B+ ^ J/rl;K*+ and B+ ^ J/rpK+ decays (7|. 

Assuming that photon reconstruction efficiencies are uncorrelated the correction factor 
is interpreted as the square of the photon reconstruction uncertainty rj))™’' = ('0^^'^) . 
This leads to following values of the corrections for the transverse energy of the photon 
more than 250 MeV. 

pV'' = (101.6 ± 1.3 (stat) ± 1.1 (syst) ± 3.3(f3)) %, 
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for 2011 data and for 2012 data is 


= (102.9 ± 0.9 (stat) ± 0.8 (syst) ± 3.4(f3)) %, 

where the third uncertainty is related to an uncertainty of the known branching fractions 
of B+ ^ J/rl;K*+ and B+ ^ J/tljK+ decays. 

7.1 Correction factor for different E^(y) 

The correction factor r]™’'’’ is also studied in four bins of the photon transverse energy. 
The correction factors can be represented as product of for hrst and second 
photons. This leads to ten correction factors for 7r°, depending on the combination of the 
bins, hit by the daughter photons of the tt*^. The system of equations is obtained and 
then solved using the method. During this procedure only statistical uncertainties are 
taken into account. 

In order to estimate the systematic uncertainties of the photon corrections, the correc¬ 
tions for the 7r° reconstruction efficiencies are varied within their systematic uncertainties. 
Each time the set of photon corrections is recalculated with the help of the method 
described above. In total 400 such experiments are performed. The 68% interval of the 
variation of the photon reconstruction efficiency correction in each bin is taken as the 
systematic uncertainty for this bin. The obtained values for r|“" are shown on Fig. and 
certain values are listed in Table [H 




Figure 4: Factors for (a) 2011 and (b) 2012 data in bins of photon transverse energy. 


7.2 Correction factor for different p^(7T®) 

The dependence of the correction factor on the neutral pion transverse momentum 
is also studied. The transverse momentum spectra are divided into a few bins and the 
analysis is repeated for each bin separately. The widths of the bins are varied to allow 
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Table 3: Factors 'n™'’'’ for 2011 and 2012 data in bins of photon transverse energy. The 
first nncertainty is the statistical nncertainty, the second is systematic. 


[GeV/c] 

ri“" [%] 2011 data 

ri™’'’’ [%] 2012 data 

250 - 350 

96.5 ±4.5 ±4.1 

102.0 ±3.9 ±3.5 

350 - 500 

96.1 ±3.2 ±3.4 

98.2 ±2.3 ±3.3 

500 - 700 

93.5 ±2.5 ±3.3 

95.4 ± 1.9 ±3.2 

> 700 

104.5 ± 1.7 ±3.6 

103.8 ±1.1 ±3.4 


for sufficient number of candidates in each bin. The obtained values for are listed in 
Table m for 2011 data and in Tables [5] for 2012 data. 

Table 4: Factors ri^o" for 2011 data in bins of neutral pion transverse momentum. The hrst 
uncertainty is the statistical uncertainty, the second is systematic, and the third is related 
to an uncertainty of the known branching fractions of —>■ J/i|) K*"*" and B+ —)■ J/t|) K+ de¬ 
cays. 


p;!, |GeV/c] ng" [%] 


0.5 

- 1.0 

93.9 

± 

00 

± 

2.5 

± 

6.1 

1.0 

- 1.25 

92.1 

± 

5.0 

± 

2.4 

± 

5.9 

1.25 

- 1.5 

94.9 

± 

4.3 

± 

2.3 

± 

6.1 

1.5 

- 2.0 

99.5 

± 

3.2 

± 

3.7 

± 

6.4 


> 2.0 

117.2 

± 

3.7 

± 

2.6 

± 

7.6 


Table 5: Factors for 2012 data in bins of neutral pion transverse momentum. The first 
uncertainty is the statistical uncertainty, the second is systematic, and the third is related 
to an uncertainty of the known branching fractions of B’*' —)• J/i|) K*"*" and B'*’ —)■ J/\|; K+ de¬ 
cays. 


p3. |GeV/cl_nS" l%l 


0.5 

- 1.0 

89.7 

± 

5.9 

± 

1.6 

± 

00 

1.0 

- 1.25 

90.6 

± 

3.4 

± 

1.5 

± 

5.9 

1.25 

- 1.5 

94.9 

± 

3.1 

± 

1.6 

± 

6.1 

1.5 

- 2.0 

104.8 

± 

2.3 

± 

1.5 

± 

6.8 


> 2.0 

116.4 

± 

2.6 

± 

1.6 

± 

7.5 
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8 Summary 

Using data corresponding to 3 fb“^, collected in 2011 and 2012 with the LHCb detector, 
the correction factors for the 7 t° and y reconstruction efficiency are obtained. Their values 
are determined through measurement of the relative yields of B"*" —>■ J/r|jK*+ (—)■ K+tt®) 
versus B+ —)■ J/t|) K"*" events. The correction factor for the photon reconstruction efficiency 
are interpreted as r|“" ~ dependency of the corrections factors on the photon 

transverse energy is studied. 
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